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Abstract. In this work, we study the Luminosity-Metallicity relation 
(LMR) and the Stellar Mass-Metallicity relation (MMR) of galactic sys- 
tems in a hierarhical clustering scenario. We performed numerical hy- 
drodynamical simulations with the chemical GADGET-2 of Scannapieco 
et al.(2005) in a ACDM universe. We found that our simulated galactic 
systems reproduce the observed local LMR and its evolution in zero point 
and slope. The simulated MMR is also in agreement with recent observa- 
tional results. From the analysis of the evolution of the MMR, we found a 
characteristic mass at ~ 10 10 ' 2 Mq /i" 1 which separates two galactic pop- 
ulations with different astrophysical properties. More massive systems 
tend to have their stars formed at z > 2 and show less evolution than 
smaller systems. Hence, this characteristic mass is determined by the 
formation of the structure in a hierarchical scenario. Our results also 
suggest the need for efficient supernova feedback. 

Resumen. En este trabajo, estudiamos la relacion Luminosidad-Metali- 
cidad (LMR) y la relacion Masa Estelar-Metalicidad (MMR) de los sis- 
temas galacticos en un modelo de agregacion jerarquica. Realizamos 
simulaciones numericas hidrodinamicas con el codigo qumhco GADGET- 
2 de Scannapieco et al. (2005) en un universo ACDM. Encontramos 
que nuestros sistemas galacticos simulados reproducen la LMR local y 
su evolucion en el punto cero y la pendiente. La MMR simulada esta 
tambien en acuerdo con resultados observacionales recientes. A partir del 
analisis de la evolucion de la MMR, hallamos una masa caracteristica en 
~ 1O 1O ' 2 M0 h , la cual separa dos poblaciones galacticas con diferentes 
propiedades astrofi'sicas. Los sistemas mas masivos tienden a formar sus 
estrellas a z > 2 y muestran menor evolucion que los sistemas pequehos. 
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Entonces, esta masa caracteristica es determinada por la formacion de 
la estructura en un universo jerarquico. Nuestros resultados sugieren 
tambien la necesidad de un importante feedback por supernovas. 



1. Introduction 

Determining the chemical composition of galaxies is of fundamental importance 
for tracing back the history of evolution of galaxies. In particular, the LMR has 
been widely studied in the Local Universe. Metallicities are tightly related with 
the luminosities of galaxies in such a way, that brighter systems have higher 
abundances (Lamareille et al. 2004). Furthermore, recent studies have also sug- 
gested that this relation extends to intermediate redshifts but displaced towards 
lower metallicities and higher luminosities (Kobulnicky et al. 2003). 
When studying galaxy evolution, stellar mass is a better parameter than lumi- 
nosity. However, because of the difficulties in obtaining stellar masses, most 
studies have used luminosity as a surrogate. Recently, though, Tremonti et al. 
(2004) have estimated the relation between metallicity and stellar mass in the 
Local Universe. The authors found a strong correlation extended over 2 dex in 
stellar mass and a factor of 10 in metallicity. 

In this work, we study the evolution of the MMR and the LMR by employ- 
ing numerical chemo-dynamical simulations which allow to describe the non- 
lineal growth of structure simultaneously with the enrichment of the interstellar 
medium in a cosmological framework. 



2. Results and discussion 

We have run numerical simulations by using the chemical GADGET-2 of Scan- 
napieco et al. (2005). A ACDM cosmological model (^=0.3, A=0.7, Q b =0.04 
and #0=100 h km s~ x Mpc -1 with h=0.7) was assumed, according to which 
galaxies formed by the hierarchical aggregation of substructures. We have anal- 
ysed two realizations of the power spectrum in a 10 Mpc h^ 1 side box, initially 
resolved with 2 x 160 3 (S160) and 2 x 80 3 (S80) particles, corresponding to mass 
resolutions of 2.71 x 10 6 M h' 1 and 2.17 x 10 7 M h' 1 for the gas phase and 
1.76 x 10 7 Mq h~ l and 1.41 x 10 8 Mq h' 1 for dark matter respectively. 
A Salpeter Initial Mass Function has been assumed with upper and lower lim- 
its of 40 Mq and 0.1 Mq, respectively. The chemical model includes metal- 
dependent radiative cooling, star formation and chemical enrichment by super- 
novae II and la (Scannapieco et al. 2005). 

Galactic objects were identified by applying an identification algorithm that 
combines the friends-of-friends technique and the contrast density criterium of 
White, Efstathiou & Frenk (1993). Dynamical and chemical properties were 
estimated at the optical radius calculated accordingly to the standard defini- 
tion as the radius which contains the 83% of the baryonic mass of the system 
(Tissera et al. 2005). Colours and magnitudes of galactic systems were calcu- 
lated by resorting to population synthesis models (see De Rossi et al. 2006 in 
preparation). 



Luminosity-Metallicity-Stellar Mass correlation in galaxies 
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Our simulations predict a linear correlation between oxygen abundance and 
luminosity which is in good agreement with the observational results. We have 
also found an evolution in the LMR in such a way that the slope increases and the 
zero point decreases with redshift consistently with the findings of Kobulnicky 
Sz Kewley (2004), among others. In particular, we decided to work with the 
I-band because it is less affected by extinction and can be more directly related 
with the underlying mass distributions. Our results indicate that at a given 
chemical abundance, galactic systems are ~ 3 dex brighter at z = 3 compared 
to z = 0, with the larger evolution at fainter magnitudes. Futhermore, we have 
encountered a mean evolution in the chemical abundances of galactic systems of 
~ 1.6 dex for brighter magnitudes and ~ 2.5 dex for faint ones, from z = 3 to 
z = 0. 

We have also analysed the MMR for simulated galactic systems, obtaining sim- 
ilar trends to those found by Tremonti et al. (2004) in the Sloan Digital Sky 
Survey (SDSS) but with a displacement of -0.25 dex in the zero point. This last 
difference may be explained taking into account that the SDSS explored only the 
central regions of galaxies which could lead to an overestimation of their metal 
content. Galactic abundances derived from simulations tend to increase with 
stellar mass which is also consistent with the observed behaviour. However, we 
obtained an excess of metals in the lower mass end which could be due to the 
absence of supernovae energy feedback in our model. 

We have determined a characteristic stellar mass at M c = 10 102 Mq h~ x where 
a change in the curvature of the MMR occurs. This characteristic mass, which 
corresponds to an oxygen abundance of ~ 8.7 dex, has been obtained by esti- 
mating where the residuals of the linear fits depart systematically from zero. It 
is important to note that this mass is similar to the characteristic mass derived 
from the SDSS by Tremonti et al. (2004) and Gallazzi et al.(2005). 
In addition, we have found that the MMR exhibits the same general patterns 
from z = 3 to z = 0, but with a displacement towards higher abundances as 
redshift decreases. The characteristic stellar mass M c remains almost unchanged 
with time and only its corresponding chemical abundance evolves by 0.05 dex in 
the same redshift range. The major departure from the local MMR occurs for 
smaller systems which increase their chemical content by ~ 0.10 dex. On the 
other hand, massive systems show less evolution with variations of ~ 0.05 from 
z = 3 to z = 0. 

We have also studied the Metallicity-optical Velocity Relation (MVR) finding a 
well defined correlation from z = 3 to z = 0. However, a higher level of evolution 
has been found in the MVR when compared to the MMR. Fast rotators show 
an enrichment of ~ 0.18 dex from z = 3 to z = while at lower metallicities 
the variations are of ~ 0.28 dex. This significant evolution of the MVR is a 
consequence of the increase of the mean density of the Universe at high redshift, 
so that at a fix stellar mass, systems need to concentrate more in order to fullfill 
the contrast density criterium and, hence, galaxies reach higher velocities. The 
larger evolution showed by the simulated MVR, when compared with the MMR, 
shows that the second is more fundamental in the sense that it is not strongly 
dependent on the cosmic epoch (see Tissera, De Rossi & Scannapieco 2005 for 
details). 
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By analysing the merger trees of systems at z = 0, we have encountered that 
the features of the MMR can be traced back in time taking into account the 
mergers and interaction history of galaxies within the hierarchical aggregation 
picture. More massive systems transform most of their gas into stars at high 
redshifts and experience important mergers. At lower redshifts, these galactic 
objects are saturated of stars, so that, while their stellar mass largely increases 
in a merger event, their metallicity remains basically unchanged. On the other 
hand, less massive systems form their stars in a more passive way or by rich-gas 
mergers leading to a more tight correlation between metallicity and stellar mass. 

3. Conclusions 

We have estimated the MMR, the LMR and the MVR correlations finding simi- 
lars trends to those detected in observations. These relations evolve with redshift 
with the major changes driven by less massive systems in consistency with the 
downsizing scenario. 

We have found an excess of metals for lower masses which could be attributed 
to the lack of an energy feedback treatment in our simulations. 
A characteristic stellar mass M c ~ 10 10,2 Mq/i _1 has been determined at which 
the MMR flattens. This mass is in good agreement with the one found by 
Tremonti et al (2004) and Gallazzi et al. (2005). 

Our results suggest that the features of the MMR are tightly related to the 
hierarchical aggregation scenario in which galaxies form. 
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